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ABSTRACT 

We construct a fully self-consistent mass model for the lens galaxy SDSS J2141 at redshift 
0.14, and use it to improve on previous studies by modelling its gravitational lensing effect, 
gas rotation curve and stellar kinematics simultaneously. We adopt a very flexible axisym- 
metric mass model constituted by a generalized NFW dark matter halo and a stellar mass 
distribution obtained by deprojecting the multi-Gaussian expansion (MGE) fit to the high- 
resolution K'-band Laser Guide Star Adaptive Optics (LGSAO) imaging data of the galaxy, 
with the (spatially constant) mass-to-light ratio as a free parameter. We model the stellar kine- 
matics by solving the anisotropic Jeans equations. We find that the inner logarithmic slope 
of the dark halo is weakly constrained, i.e. y = 0.82^q^4, and consistent with an unmodified 
NFW profile; we can conclude, however, that steep profiles (y ^ 1.5) are disfavoured (< 14% 
posterior probability). We marginalize over this parameter to infer the galaxy to have (i) a 
dark matter fraction within 2.2 disk radii of 0.28^q|q, independent of the galaxy stellar popu- 
lation, implying a maximal disk for SDSS J2141; (ii) an apparently uncontracted dark matter 
halo, with concentration c_2 = ^.T^^'s virial velocity v^ir = 242^3^ kms"\ consistent 
with ACDM predictions; (iii) a slightly oblate halo (qh = 0.75^^^^), consistent with predic- 
tions from bary on-affected models. Comparing the tightly constrained gravitational stellar 
mass inferred from the combined analysis (log^oM^/M© = 11. 12^^^^) with that inferred 
from stellar populations modelling of the galaxies colours, and accounting for an expected 
cold gas fraction of 20 ± 10 per cent, we determine a preference for a Chabrier IMF over 
Salpeter IMF by a Bayes factor of 5.7 (corresponding to substantial evidence). We infer a 
value Pz = I- o'I/o'r = 0-43 ^^^^ for the orbital anisotropy parameter in the meridional plane, 
in agreement with most studies of local disk galaxies, and ruling out at 99 per cent confidence 
level that the dynamics of this system can be described by a two-integral distribution function. 

Key words: gravitational lensing: strong — galaxies: fundamental parameters — galaxies: 
haloes — galaxies: kinematics and dynamics — galaxies: spiral — galaxies: structure 
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1 INTRODUCTION 

Measuring the relative contribution of luminous and dark matter in 
spiral galaxies is essential to understand their internal structure and 
therefore constrain the physical process es that drive their formation 
and evolution (e.g. IPutton et alJEoilah . Traditionally this is done 
by means of detailed s tellar and gas kinema tics and stellar popu- 
lation diagnostics fe.g. lBershadv et al1l201lh . However, often one 
needs additional assumptio ns about the relat ive contribution of the 
stars and dark matter (e.g. Ivan Albada & Sancisi 1986). or about 
the stellar initial mass function (e.g. lBell & de Jondl200lh . 

The combination of strong gravitational lensing and galaxy 
kinematics is a powerful tool for constraining the mass distri- 
bution and the dynamical structure of galaxies beyond the lo- 
cal Universe, since this approach makes it possible to overcome 
many of the difficulties associated with the traditional techniques, 
whi ch are severely lim ited when app lied to dis tant objects (see 
e.g. lTreu&KooDmanslt2002, 2004; Barnabe & Koopmans 20071 ; 



Jiang & Kochanekl l2007iriGrillo et al.l boOSi ; lAuger et al.ri2010al ; 
Koopmans et al.l l2009l ; iBarnabe et all l20ld l201lh . In particular. 



gravitational lenses in which the deflector contains a high incli- 
nation disk provide extra (geometrical) information to help disen- 



the 3D shape of the dark matter halo (e.g., Keeton & Kochanekl 


1998L Koopmans et a 


.1 19981 iMal 


er et all 200d iTrott et al.ll2010l. 


Suvu etal.ll201ll and 


Dutton et al. 


l2011bl. Paper II of this series) 



Because these measured masses are gravitational, they can be com- 
pared with the stellar mass from stellar population synthesis (SPS) 
models and so used to constrain the form of the stellar initial mass 
function (IMF), and the response of dark matter haloes to galaxy 
formation. 

Until recently only a small number of gravita tional lenses 
with high-in clination disks were known. The SLACS (iBolton et al.l 
I2OO6LI2OO8I) and SWELLS (iTreu et al.ll20Tll . Paper I) surveys have 
significantly increased the number of known gravitational lenses in 
which the deflector contains a high-inclination disk, including sev- 
eral disk-dominated systems. 

One of the most promising spiral lens systems for a joint lens- 
ing and dynamics analysis is SDSS J2141-0001 (hereafter simply 
referred to as SDSS J2141 for brevity), at redshift ziens = 0.1380, 
which is a disk-dominated galaxy (it has a disk K'-band light frac- 
tion of ^ 80%) at high inclination (/ ^ 78°). In addition to the dis- 
covery data from the SLACS survey, a wealth of imaging and kine- 
matic data are available from the SWELLS project (Paper I and II). 
A joint strong lensing and gas kinematics (rotation curve) analysis 
of SDSS J2141 conducted in Paper II yielded a gravitational stel- 
lar mass of logio(M^/Mo) = 10.99+Q 25 (consistent with that from a 
stellar population analysis assuming a Chabrier IMF), a dark matter 
fraction at 2.2 disk scale lengths of /dm = O-^S^q^^, and a dark mat- 
ter halo flattening of = 0-91^q However, in that work, simple 
phenomenological ("Chameleon") models for all three mass com- 
ponents, i.e. the dark matter halo and the stellar bulge and disk, 
were assumed. Moreover, only a fraction of the available kinematic 
data for the lens galaxy was used: the stellar velocity dispersion and 
rotation curve were not considered. Indeed, the velocity dispersion 
could not be predicted self-consistently within the assumed model. 

In this paper we improve on the Paper II analysis of 
SDSS J2141 in several important ways. The main improvement is 
the inclusion of stellar kinematics data, which provides a mass con- 
straint at smaller radii than obtained from lensing or gas kinemat- 
ics. It is well known that disk galaxies are usually characterized 
by a velocity dispersion ellipsoid flattened along the vertical direc- 



tion: therefore, in order to provide an accurate description of the 
data set, we model the stellar kinematics by means of anisotropic 
Jeans equations, which allow us to properly take into account (and 
recover) the anisotropy ratio parameter p.. In addition to this, we 
use more flexible and general models for both the stellar and the 
dark matter density profiles. Specifically, we obtain the stellar mass 
density profile from the deprojection of the observed luminous dis- 
tribution (fitted with the state-of-the-ar t method of multi-Gaussian 
expansion, MGE, see ICappellarill2002l) . rather than the sum of two 
"Chameleon" profiles, wh ich were used to approximate a Sersic 
profile bulge (ISersid ll968) and an exponential disk. Finally, here 
we model the dark matter halo with an ellipsoidal generalized NEW 
profile (inner logarithmic slope -7, outer slope -3) rather than the 
non-singular isothermal ellipsoid (inner slope 0, outer slope -2) 
with a fixed density profile in the inner regions used in the previous 
analysis. 

The resulting model is both self-consistent and, in the case of 
the dark matter halo, physically-motivated, and allows us to attempt 
to fit all the data we have for SDSS J2141 simultaneously. We use 
it to answer the following questions about SDSS J2141: How much 
does dark matter contribute to the total mass of this disk galaxy, in 
particular in the inner regions? What is the concentration and in- 
ner profile slope of its dark matter halo? What is its halo's shape? 
When calibrated via its stellar mass distribution's gravitational ef- 
fects, what galaxy- averaged IMF do we infer from a stellar popu- 
lation synthesis analysis of its optical and near infra-red colours? 
What is the vertical-to-radial anisotropy of its velocity dispersion 
ellipsoid? 

This paper is structured as follows. We first describe our obser- 
vational data (imaging for the lens modelling, spectroscopy for the 
stellar and gas kinematics) in Section |2] We then outline our mass 
model for SDSS J2141 in Section[3] giving the functional forms we 
use to describe its stellar and dark matter distributions. Then, in 
Section |4] and Section |5] we show how our model predicts both the 
lensing and kinematic data in a self-consistent way, and in Section[6] 
we review the probability theory behind the actual inference proce- 
dure we follow. In Section |7] we present and discuss the results of 
our analysis, and in Section [8] we draw conclusions, providing an 
answer to the questions posed above. 

Throughout this work, we assume a flat ACDM cosmology 
with present day matter density, Qm = 0.3, and Hubble parameter, 
Hq = 70kms-iMpc. 



2 OBSERVATIONS 

In this Section we briefly recall the data set available for this study. 
A more detailed description of the data is given in Paper II. 



2.1 Imaging data 

The imaging data consists of a high spatial resolution (FWHM 
^ 0.15 arcsec) K'-band image taken with adaptive optics on the 
Keck II telescope. The galaxy- subtracted image (see Paper II) is 
used for the strong lensing analysis, while the light profile of the 
galaxy is fitted with a set of elliptical Gaussians which are de- 
convolved and deprojected to provide a 3D model of the stellar 
mass (up to the normalization), as detailed in Section [3Jl The lens- 
subtracted image used as data set for the lensing analysis is shown 
in the upper-right panel of Figure [T] Multi-band HST photometry 
is also available and used to determine stellar mass as discussed in 
Papers I and II. 
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Figure 1. Lensed image reconstruction obtained from the maximum a pos- 
teriori model. From the top left-hand to bottom right-hand panel: recon- 
structed source; HST/ACS data showing the observed lensed image after 
subtraction of the lens galaxy; lensed image reconstruction corresponding 
to the source in the first panel; residuals. 



2.2 Kinematic data 

The second set of data that we will use to constrain our mass model 
is the rotation and velocity dispersion profiles derived from opti- 
cal emission- and absorption-line spectroscopy. A major axis long- 
slit spectrum of SDSS J2141 was obtained with the DEep Imaging 
Multi-Object Spectrograph (DEIMOS) on the Keck II 10-m tele- 
scope. 

We used the 1200 line grating (corresponding to a pixel scale 
of 0.32 A) with a 1 arcsec width slit resulting in a spectral resolution 
of ~ 1.9 A. The wavelength range was 5200 - 7800 A, covering 
several prominent emission and absorption lines. At the wavelength 
of Mgb the velocity resolution is cr^^s - 41 kms~^ , and for Ha it 
is (Tres - 32 km s~^ . We took three exposures of 1200s in excellent 
seeing conditions of 0.60 arcsec. The spectra were reduced using 
routines developed by D. Kelson (lKelsonll2003h . 

Kinematic parameters were measured from one-dimensional 
spectra extracted along the slit with a spatial sampling of ^ 0.59 
arcsec (5 DEIMOS pixels), corresponding to one data point per 
seeing FWHM. The rotation and velocity dispersion profile of the 
stars were obtained by fitting a region including the Mgb [5177 A] 
and Fell [5270 A] lines with a set of stellar templates. 

The rotation curve of the ionized gas was measured by fitting 
Gaussians to the Ha line [6563 A], and is shown in the upper panel 
of Figure [2 (data points with error bars). Outside of the inner ~ 2 
arcsec the velocity dispersion of the Ha line was equal to the instru- 
mental resolution, indicating the ionized gas disk is dynamically 
cold. 

In our dynamical model (see Section [5]), we assume that the 
ionized gas traces the circular velocity of the galaxy (i.e., there is no 
pressure support). For the stellar kinematics our model implicitly 
includes rotation and dispersion, although neither of these parame- 
ters are fitted to directly. Instead, our model predicts the projected 
second velocity moment, which is fitted to the root mean square 
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Figure 2. Observed gas and stellar kinematics compared to the predictions 
of the maximum a posteriori model. The upper panel shows the galaxy Ha 
rotation curve (blue data points): the red line represents the intrinsic model 
circular velocity, while the black line gives the predicted observable, i.e. 
the model circular velocity after the beam- smearing, finite slit width and 
inclination effects are taken into account. The light blue data points are not 
used to constrain the model (see text). The lower panel shows the model 
projected second velocity moment 112 (black line) compared to the corre- 
sponding observational quantity f;rms (blue data points). The stellar kine- 
matic data sets |i;rotl (stellar rotation curve, green) and cr (stellar velocity 
dispersion profile, red) are also shown for reference. See Section |5] for a 
rigorous definition of these quantities. 



velocity of the stars, v^^{R) - Jv^^^(R) + o-^(R) f see ICappellaril 
l2008h . The lower panel of Figure|2]shows the observables v^ot and cr 
(as green and red data points, respectively) as well as the root mean 
square velocity (blue data points). 

In our modelling, as described in Paper II, we conservatively 
exclude from the fit the data points of the gas rotation curve that 
are (i) within the inner 2 arcsec (due to uncertainties and likely 
asymmetries in the Ha distribution in this region) and (ii) beyond 
3.5 arcsec on the west side of the rotation curve (where there is an 
asymmetry caused by the presence of the warp). These excluded 
points are shown in light blue in Figure |2] 



3 THE GALAXY MASS MODEL 

In order to perform a self-consistent analysis of the mass structure 
of SDSSJ2141 we need to combine the constraints derived from 
both the lensing and kinematics data sets. The most general and 
straightforward way to proceed is simply to adopt for the galaxy a 
plausible total mass density distribution ptot(-^, ^), where x denotes 
the spatial coordinates and // is a set of parameters characterizing 
the density profile, and use it to model simultaneously the various 
sets of observables. The main challenge with this approach lies in 
choosing a mass distribution that is realistic and flexible enough 
to reproduce the data, but at the same time simple enough that the 
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exploration of the 1/ parameter space remains computationally fea- 
sible. 

In keeping with local studies of both disk and e lliptical galax- 
ies fe.g. lWeiimans et al.|[2008L IWilHams et alJl2009l) we model the 
mass distribution using two components: (i) a luminous mass com- 
ponent whose detailed profile is obtained by deprojecting the ob- 
served surface brightness distribution and (ii) a generalized NFW 
dark matter halo whose profile is motivated by cosmological simu- 
lations. 

Throughout the remainder of this work we will make reference 
to the two following right-handed coordinate systems: (i) a cylin- 
drical coordinate system {R, ^, z) with the z-axis directed along the 
galaxy rotation axis; (ii) a Cartesian coordinate system {x\y\z'), 
where the z'-axis is directed along the line-of-sight and {x\ y') de- 
notes the plane of the sky, with the x'-axis aligned along the galaxy 
projected major axis. In both cases the origin of the axes is located 
in the center of the galaxy (which is assumed to be an axially sym- 
metric system). We use the first reference system to write the in- 
trinsic galactic quantities (e.g. density and potential) and the sec- 
ond one to express the projected quantities (e.g. surface brightness 
and deflection angle). We denote with / the galaxy inclination, i.e. 
the angle comprised between the z- and z'-axes (so that / = 90° for 
a system observed edge-on). 



3.1 Luminous mass distribution 

An ideal model for the luminous mass distribution should be 
flexible and realistic (in particular, it should be able to reproduce 
the observed surface brightness distributions when projected 
along the line-of-sight), and analytically simple, so that the 
corresponding gravitational potential is easy to calculate. This 
can be achieved by makin g use of the MGE method, a technique 
originally pioneered by iBendinelhl (Il99ll) and subsequently 
genera lized and develop ed bv [Monnet Bacon. & Eniselleml 

(Il992h. [Emsellem. Monnet & BaconI ( Il994l) . 

lEmsellem. Deionghe. & BaconI (Il999l) and Cappellaril (l2002h , 
whose formalism we follow here. In order to minimize dust 



obscuration and map as closely as possible the stellar distribution, 
we apply the MGE decomposition to the high-resolution K'-band 
image of the SDSS J2141 surface brightness distribution. 

The observed galaxy surface brightness l.{x' ,y') is 
parametrized as a sum of N two-dimensional, concentric, 
elliptically- symmetric Gaussian components g^ix'.y'), each with 
luminosity L^: 



I^(x\y') = ^Lkgk(x\y'), 

k 

where each Gaussian function 
1 



gk(x\y') 



exp 



(1) 



(2) 



is characterized by the widths cr^ and q'^ak along the x'- and ^'-axis 
respectively, and q'^ is the projected axial ratio of the ^-th compo- 
nent. The total stellar luminosity of the system is simply given by 

In general, even assuming — as we do — that the galaxy in- 
clination angle / is known, the deprojection of the observed light 
distribution of an axisymmetric galaxy is an intrinsically degen- 
erate problem unless the system is seen edge-on (lRvbickilll987h . 
The solution, however, can become unique when a model is speci- 
fied. In the case of the MGE parametrization, the deprojected three- 



dimensional luminosity density distribution has the simple expres- 
sion 



{Infl^o-lq, 



exp 



1 









(3) 



which is still a sum of Gaussian functions with intrinsic axial ratios 
given by 

^ - cos^ / 

-T—- (4) 



^k 



Since the stellar component of galaxies is oblate or spherical, most 
(if not all) Gaussians will turn out to have ^ qk ^ 1- In order 
for the axial ratios of the 3D Gaussian components to be physical, 
one must enforce the constraint that the projected axial ratios q[ are 
rounder than |cos /| when fitting the profile of Eq. ([U to the observed 
surface brightness distribution. 

The luminosity density in Eq. ^ can be straightforwardly 
converted into a mass density by multiplying each term by the stel- 
lar mass-to-light ratio T^, so that the mass of each Gaussian is given 
by = T^Lyt. However, the single Gaussian elements are simply a 
mathematically convenient way to describe the light profile and do 
not have a direct physical meaning individually. Therefore, since 
there is little interest in studying them one by one, in this work we 
assign the same global stellar mass-to-light ratio to all the lu- 
minous components. This simplifying assumption is equivalent to 
assuming that the bulge and disk components have the same stellar 
mass-to-light ratio. Note that the choice of the K'-band image as 
trace of stellar light should minimize variations in mass-to-light ra- 
tio. Future work with higher resolution data should explore further 
the limitations of this assumption. 

An additional advantage of the MGE approach is that we avoid 
dealing with the difficult and somewhat degenerate problem of de- 
composing the lig ht profile into the separate dis k and bulge con- 
tributions (see e.g. Ivan der Kruit & Searlelll98lh since we have a 
model that can fit very accurately the whole light distribution at 
once. 

The density distribution of Eq. (|3]l is a sum of components 
stratified on homoeoidal surfaces, hence the corresponding gravi- 
tational potential can be deriv ed using the cla ssic Cha ndrasekharl 
(Il969h formula, obtaining (see lEmsellem et al.|[T994l) 



0(/?,z) 



(5) 



where G is the gravitational constant and the dimensionless func- 
tion 



Jo 



dr 



: exp 



(6) 



with ri\ = \ - q\, can be evaluated with a single numerical inte- 
gral. The density distribution (Eq. [3]) and its potential (Eq. [5]) are 
remarkably simple for such a flexible mass model. Even better, the 
corresponding kinematic quantities, obtained by solving the Jeans 
equations, also have relatively straigh tforward expressi ons that do 
not involve any special functions (see ICappellarill2008l for a rigor- 
ous derivation of the velocity moments). 

3.2 Dark matter halo 

Cold dark matter simulations are known to produce halos with, on 
average, universal mass density profiles dNavarro, Frenk. & White! 
[1997|, NFW) that are well fitted by a broken power-law functional 
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form with an inner logarithmic slope y - \ and a slope 7 = 3 in 
the outer regions, i.e. at radial distances much larger than the scale 
radius r^. The situation, however, becomes far more complex when 
the baryons are added to the picture and — although the detailed 
mechanisms are not yet fully understood — it is widely accepted 
that the involved processes can have the effect of modifying the in- 
ner s l ope of the dark ha l o density profile (e .g. iBlumenthal et al.' 
19861: iDekel et all l2003l : iGnedin et all l2004l : iNipoti et al. 2004,; 
Abadi et al.ll201Qh . Therefore, in order to account for a dark matter 



halo that can be either steeper or shallower than a NFW in the in- 
ner regions, we adopt a n axisymmetric generalized NF W (gNFW) 
density distribution (see lZha ol ll996l : IWvithe et al.ll200lh : 



Pdm(w) 



Pa 



(m/ny (1 + m/rs)3-r ' 



(7) 



Here, pcht is the critical density of the universe at the redshift of the 
object in question, and m denotes the elliptical radius, i.e. 



(8) 



where indicates the three-dimensional axial ratio of the pro- 
file (the halo is oblate for < 1, spherical for q\^- I, prolate for 
> !)• Note that we take the halo to be aligned with the stellar 
mass distribution, as in Paper II. 

In order to enable an easier comparison of the scale radii be- 
tween profiles having different values of 7, it is useful to introduce 
the quantity r_2 = (2 - y)r^, which corresponds to the radius at 
which the logarithmic density slope of the profile is -2. Clearly, 
r_2 = Ts only in the case of the regular NFW profile. Another use- 
ful scale length is the "virial" radius ryir, defined as the spherical 
radius within which the average density is equal to 200 pcrk- The 
concentration parameter of the halo is usually expressed as the ra- 
tio c - Tyir/r^', an alternative definition, adopting the radius r_2, is 

The characteristic halo density 6c that sets the normalization 
of Pd m in the center is then expressed (following, e.g. JPutton et al.l 
I2005li1 as a function of both the concentration and the slope: 



6c ^ 



200 



3 f(c,7,l)' 
where we have defined the function 



f(c,7,^h) ^ 



Jo ^1 -(1-^2)^2/, 



:dT. 



(9) 



(10) 



The mass distribution given by Eq. is completely specified 
when the four independent parameters 7, q^, and c are given. 
In this work, we choose to re-parametrize the halo using the virial 
velocity Uyir, i.e. the circular velocity at the virial radius, in place 
of the scale radius, since Uvir has a very intuitive physical interpre- 
tation and facilitates the comparison with the oretical work, where 
this quantity is freq uently employed (see e.g. iMaccio et al.ll2008l : 
button etal.ll2011ah . If the velocity is expressed in km s ' and the 
radii in kpc, then one can show (cf. .Dutton et alji2005i^ that Oyir is 
related to the virial radius by the formula 



(11) 



^ Note that there is a typographical error in Eq. (7) of IPutton et all l2005h : 
inside the integral the numerator should read y^~^[l + (2 - a)y]^~^. 



where h = ///lOOkms"^ Mpc~^ and H denotes the value of the 
Hubble constant at the redshift of the object. 

It is convenient to perform an MGE of the axisymmetric 
gNFW profile in order to simplify considerably both the calcula- 
tion of t he lensing angle and the solving of the Jeans equations 
(cf., e.g. JWilliams etal.ll2009l . where a MGE decomposition of the 
NFW halo is performed). We find that around 8 Gaussian com- 
ponents are typically enough to provide an excellent fit to both the 
Pdm distribution and the lensing deflection field (typically within 1- 
3 per cent), ensuring that the adoption of this approximation does 
not change our inferences. In this case, the total potential is still ob- 
tained from Eq. ^ by extending the sum to include also the A^dm 
Gaussian elements that describe the dark halo component. 



4 MODELLING THE GRAVITATIONAL LENSING 

Given the observed surface brightness distribution of the lensed 
images and a mass model for the deflector, we recover the un- 
lensed surface brightness distribution of the background object (the 
"source" object) by making use of the pixellated source reconstruc- 
tion method, which takes i nto account the effec t s of PSF blurring 
and regularization (see e.g. Warren & Dvel l2003l : lKooDman"s 2005: 
ISuvu et aP I2OO6I : iBrewer & Lewis 2006). Our implementa t ion of 



Barnabe & KoopmansI (l2007h 



this method is described in detail in|] 
and is included in the CAULDRON code that has been employed in 
the combined lensing and dynamics analysis of the SLAGS early- 
type gala xies for which two-d i mensional kinematic maps a r e avail- 
able (see ICzoske et al.ll2008L l2012l : iBarnabe et al.ll2009al . l2009bl . 

Imil). 

This approach consists in casting back, pixel by pixel, through 
the lensing equation, the lensed image grid onto the source image 
grid. The results of this procedure are encoded in the lensing oper- 
ator A, which allows one to express the mapping of the background 
source s into the lensed image as a linear problem, i.e. As = d. 
This set of linear equations is then solved for s by means of very 
efficient standard techniques. 

All that is needed to calculate the lensing operator is 
the deflection angle or, which is obtained from the sur- 
face mass density distribution of the le ns galaxy (see e.g. 
ISchneider. Kochanek. & WambsganssI l2006h . Therefore, the ma- 
trix A depends both on the physical parameters rj characterizing 
the density profile and on the geometry of the system, i.e. the in- 
clination / and the angular diameter distances between the observer 
and the source (D^), between the observer and the lens (D^), and 
between the lens and the source (D^s)- 

For many three-dimensional density profiles of astrophys- 
ical interest the deflection angle is very cumb ersome to com- 
pute (cf., e.g., the catalogue of iKeetonI l200lh . This has con- 
tributed to the widespread adoption of those few profiles, such 
as the isothermal ell ipsoid, for which analytical expressions 
of a are available (see Kormann, Schneider. & Bartelmann|[l994l : 



iKeeton & KochanePll998h . Remarkably, the lensing deflection an- 
gle corresponding to the density distribution of Eq. ^ is very 
straightforward to calculate, involving a single quadrature and no 
special functions: 



7r/)^Scrit 



Jo 



TdrV 



X exp 



2 \ 1 



(12) 
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f 

Jo 



rdrV 



xexp|-y 



3/2 



(13) 



where both the deflection angle and the widths cr^ are expressed in 
radians; x' = x' IcTk and y' = y' IcTk are the sky coordinate (normal- 
ized to (Tk) with respect to the lens center. As before, t/^ = (1 - ^^). 
The critical surface density 



(14) 



is th e characteri stic surface density used in gravitational lensing 
(e.g.,Ei3[20i3). 



5 MODELLING THE KINEMATICS 

5.1 Predicting the observed second velocity moments with 
the anisotropic Jeans equations 

Let us consider a steady- state axially symmetric stellar system 
characterized by a distribution function (DF) /(x, v), where the po- 
sitions X and the velocities v are the phase-space coordinates, and 
subject to the influence of a total gravitational potential Otot(^,z). 
While the typical observational data sets do not allow one, in gen- 
eral, to recover the full six-dimensional DF, it is possible to gain 
valuable information on the global dynamical structure of the sys- 
tem by no ting that its velocity mome nts must satisfy the two Jeans 
equations teinnev & Tremainell2008l) 



d(pvl) 1 d(RpU^ 
dz R OR 



= -P 



d(pvl) 



dR 



dz 



dz 
dR 



+ p- 



R 



(15) 
(16) 



Here, p{R, z) = f fd^v denotes the three-dimensional density 
distribution of the stellar system, and the bar indicates a phase- 
space average of the quantity of interest, i.e. 



(17) 



The system is not required to be self-gravitating and therefore 
in Eqs iTSl and ([T6l l p might well be the density distribution ptr 
of a tracer stellar component described by a DF /tr and subject to 
an external potential. Moreover, if other collisionless components 
(each one defined by its own DF) are present, each one will obey 
its own set of Jeans equations within the same total potential Otot- 

In our dynamical analysis of late-type galaxies, we adopt an 
axisymmetric total potential O^t = + ^dm, where the two com- 
ponents represent the potentials of the luminous distribution and 
dark matter halo, respectively. We then write down and solve the 
Jeans equations, using the stellar density distribution p^ associated 
to the corresponding potential via the Poisson equation, in order 
to obtain the intrinsic velocity moments H These are then projected 
along the line-of- sight and — after taking into account the effect 



of instrumental PSF and aperture integration — compared with the 
corresponding observational quantities. 

Of course, even when the potential and the density distribu- 
tions are given, the two Equations {15]) and ( [T6l ) still depend on the 
four unknown functions v\,v^^, vj and vrv., and therefore additional 
assumptions are needed in order to determine a unique solution for 
the Jeans equations. This is usually achieved by assigning the ori- 
entation and the shape of the intersection of the velocity dispersion 
ellipsoid with the meridional plane (R, z) at each point. 

Observations of the Milky Way and of nearby disk galax- 
ies show that the velocity dispersion ellipsoid is more flat- 
tened in the verti cal direction than in the radial onj^ i.e. 

< (se e e.g. IWielenI 1977. iGerssen. Kuiiken. & Memfieldl 

19971 bood Ivan der Kruit & de GriisI Il999 ' 



Shapiro. Gerssen. & van der Marel' '2003'. 'Ciardul lo et al.] 12004 
Noordermeer. Merrifield, & Aragon- Salamanca 2008 



ICappell^ (l2008h introduced a simple and effective way (re- 
ferred to as anisotropic Jeans models) to provide a closure for the 
Jeans equations that manages to reproduce this important feature. 
The two assumptions of this phenomenological model are: (i) the 
velocity dispersion ellipsoid is aligned with the cylindrical coor- 
dinate system (so that the mixed terms vrv^ are everywhere zero) 
and (ii) the anisotropy in the meridional plane is constant, i.e. 

= bvl, with the anisotropy parameter b ^ 0. The meridional 
plane anisotropy is usually expressed in the literature using the pa- 
rameter fiz, such that 



1 - = 



(18) 



In real galaxies, the shape and the orientation of the velocity 
dispersion ellipsoid are in general a non-trivial function of the po- 
sition on the meridional plane. However, the assumption of cylin- 
drical alignment is quite accurate for fast-rotating galaxies and disk 
systems in general, in particular along the minor axis and, more cru- 
cially, in the vicinity of the equatorial plane, where the density is 
highest. In fact, Jeans models constructed with this simple prescrip- 
tion for the anisotropy have been shown to reproduce remarkably 
well the observed kinematic moments of fast rotators a nd spirals 
dScott et al.P 2009. WiUiam s. Bureau. & Ca ppellari 2009). 

With these assumptions, the Jeans Equations ([T5] | and ([T6l ) be- 



dipyj) 

dz 
dipv{) 



+ p- 



(19) 



(20) 



dR ^ dR ^ R 

From the equations above, and imposing the intuitive con- 
straint that the vertical pressure pvj - for z oo, one obtains 
the following expressions for the intrinsic second velocity moments 
along the coordinate directions: 



^ 1 r d^,,. , 



dR 



dR 



(21) 



(22) 



^ Because of the collisionless nature of dark matter, one could write an 
analogous set of equations also for the halo component. However, since 
the corresponding velocity moments cannot be observed, this would be of 
no use in the present context of comparing the model predictions with the 
observed data sets. 



^ If the assumption of a steady- state axisymmetric system holds, this im- 
plies that the disk DF also depends on a third, non-classical, integral of 
motion I^, in addition to the two classical integrals, namely the energy E 
and the angular momentum along the rotation axis. 



© 0000 RAS, MNRAS 000, 000-000 



The SWELLS survey. IV 7 



These intrinsic quantities are then integrated along the Hne of sight 
to obtain the projected second velocity moment v^^^ (whose square 
root is usually indicated as 1^2) which can be directly compared to 
the stellar kinematics observables. The projected velocit y moments 
for the case of MGE parametrization are calculated in ICappellaril 
(l2008h . The observational counterpart of the model quantity /U2 is 

the root mean square velocity Urms = -yj^mt + where v^ot and 
cr indicate the line of sight projected stellar rotation velocity and 
velocity dispersion, respectively. 

We recall that, given a potential and a density distribution, 
the Jeans equations impose a condition for the equilibrium on 
the second velocity moments, but they provide no prescription 
on how to separate these moments into the contributions of ran- 
dom and ordered motions. Since jio net radial or vertical motions 
are considered, and thus cri = vl and cr? = u?, here this issue 
would only be relevant for the splitting of the azimuthal compo- 
nent into the streaming motion and the velocity dispersion cr^, 
i.e. = v^^ + cr^, whi ch is u s ually tackled by adopting ad hoc as 



sumptions such as the ISatohl (1 19801) decomposition. In this work, 
however, we avoid making any additional assumptions in order to 
model separately, and we only model the second velocity mo- 
ments as described above. 



5.2 Predicting the observed gas circular velocity 

In order to model the rotation curve of the Ha gas, we calculate the 
circular velocity profile, vdR), of a test particle of negligible mass 
in a circular orbit in the equatorial plane of the galaxy. 

The circular velocity, as it is clear from its definition, i.e. 
vl(R) = R (50tot/^^)lz=o. is uniquely determined by the total gravi- 
tational potential of the galaxy. In general, Vc differs from the stellar 
rotation velocity (often referred to as streaming motion) which 
is usually lower due to the effect of t he stellar velocity dispersio n, 
which acts as a pressure term (see e.g. lBinnev & TremaindbOOSh . 

The (squared) circular velocity profile that corresponds to the 
MGE mass model described in Section[3]is readily calculated from 
Eqs Is])-© and has the following expression: 



VI(R) : 



47tGR- 



qkPok 



4 



: exp 



(23) 



where again tjI = (I - q\) and the constant po^t = pyt(0,0) is the 
central value of the mass density distribution of the ^-th Gaussian 
element. 

In order to compare the predicted rotation curve vdR) of 
Eq. (I23] ) with the observations, we also take into account the com- 
bined effects of inclination, PSF blurring and finite slit width, col- 
lectively referred to as beam- smearing. Since the exact distribution 
of the Hof gaseous component is not known from the observations, 
we approximate it using the available ^'-band light profile instead, 
which is more accurate than using an exponential disk model. Ad- 
ditionally, in keeping with Paper II, we have cautionarily excluded 
from the fit two regions of the rotation curve: (i) the inner 2 arcsec, 
a region where the Ha distribution is probably asymmetric, due to 
the effect of extinction and (ii) the outermost three points of the 
west side of the rotation curve, which are affected by a spurious 
decrease of the velocity caused by the warp. 



6 BAYESIAN INFERENCE AND UNCERTAINTIES 

In order to derive rigorous constraints on the parameters that char- 
acterize the adopted model, we conduct our analysis within the 
stan dard framework of Bay esian statistics (see, e.g., Mac Kay 20031 
and ISivia & Skillh^l2006l for an extensive treatment of this sub- 
ject). 

Let us denote the combined data sets as d and the considered 
hypothesis as l-i{9). In our case, for instance, includes the model 
that we have adopted to describe the mass distribution and dynam- 
ics of the galaxy under study (Sections [3}l5]), and also all the as- 
sumptions we make about the uncertainties on the data, instrument 
response functions and any prior knowledge of the situation we 
might want to include. The non-linear parameters 6 may include, 
in general, not only the physical parameters 1/ defining the total 
mass density distribution, but also the parameters that characterize 
the dynamics (e.g. anisotropy) and the geometry (e.g. inclination) 
of the system. 

From B ayes' theorem, the posterior probability distribution 
function (PDF) for the set of parameters 6 is given by 



Vv{9\d,'H) = 



vv{d\e,9{)Vv{e\n-i) 

Vx{d\l-i) ' 



(24) 



where Pr {d \ 6, IH) is the likelihood, Pr {6 \ 7Y) is the prior, and 
Pr {d I TY), i.e. the factor required to normalize the posterior over 9, 
is known as the Bayesian evidence, which is used in comparing 
different model forms. When modelling the lensing and kinematic 
data we do not keep track of the value of the evidence, but do make 
use of it in S ection VT^ below. 

The set of parameters ^map for which the posterior is max- 
imized identifies the maximum a posteriori (MAP) model. The 
MAP model can be interpreted as a "best model" of sorts, in the 
sense that it represents the combination of parameters that is found 
to best reproduce the data given our assumptions. We adopt it 
as our reference model for the times when we need to show our 
best estimates of the predicted observables (lensed image, rotation 
curve, velocity moments) and the reconstructed background lensed 
source. 

The primary quantities of interest are the marginalized pos- 
terior PDFs for individual parameters 6i obtained by integrating 
the joint posterior PDF over all the other parameters. These inte- 
grals can be performed most readily if we characterise the joint 
posterior by a set of sample parameter values drawn from it. The 
marginalised distributions are then readily approximated by his- 
tograms of these samples. When a more compact representation 
is required, we quote parameter constraints as the median values 
of these one-dimensional histograms, ^med, and quantify our uncer- 
tainty with their 68% credible intervals (CIs, calculated by taking 
the 16th and 84th percentiles). 

The model that we employ for the analysis of SDSSJ2141 
has six free parameters, i.e. parameters with uninformative priors 
which are allowed to vary and for which the posterior exploration 
is performed. These are: the virial velocity Uvir, the inner loga- 
rithmic slope 7, the concentration c_2 and the three-dimensional 
axial ratio which describe the gNFW dark matter halo (Sec- 
tion [3]2])^ the global stellar mass-to-light ratio of the luminous 
component (Section lTTI ), that is more readily interpreted when ex- 
pressed in terms of the total stellar mass = T^Ltot; and the 
meridional plane anisotropy parameter b (Section ISTI l. In analogy 
with Paper II, we adopt a broad Gaussian prior for Uyir centered on 
255 km s~^ with a width of 45 km s~^ This corresponds to the 
prior adopted in Paper II for the virial velocity of their non- singular 
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Figure 3. Marginalized two-dimensional (contour plots) and one-dimensional (histograms) posterior PDFs for the galaxy model parameters using constraints 
from the gravitational lensing data set only. The three contours indicate the regions containing, respectively, 68%, 95% and 99.7% of the probability. 



isothermal halo, and is equivalent to assuming that the scale radius 
is not so large that the virial velocity is dramatically larger than the 
observed rotation velocity. We also adopt, as in Paper II, a lognor- 
mal prior centered on = 1 (spherical) for the axial ratio, which 
allows for both oblate and prolate haloes. We let c_2 vary on a wide 
uniform prior from to 50 to represent our ignorance of the halo 
concentration. The inner logarithmic slope y is allowed to vary be- 
tween (flat core) and 2 (isothermal). We note that Uvir is allowed 
to go to zero, so that we are also including in our analysis the case 
of a disk galaxy with no dark matter halo, fully described by a self- 
gravitating stellar mass distribution. By letting the total stellar mass 



vary (with uniform prior) between and 5 X 10^^ M© we allow for a 
wide range of contributions of the luminous components to the total 
mass, including the limiting case in which the galaxy is fully dark 
matter dominated and the stars are only a tracer with negligible 
mass. Finally, the anisotropy parameter can vary uniformly from 
^ = 0, indicating a velocity dispersion ellipsoid without any radial 
component, to Z? = 5, for which the velocity dispersion ellipsoid 
is very elongated along the radial direction: this interval is wide 
enough to include all the values of meridional plane anisotropy ob- 
served in real disk galaxies (see Section 1731 1. The model parame- 
ters, together with the adopted priors, are summarized in Table [T] A 
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Figure 4. Marginalized two-dimensional (contour plots) and one-dimensional (histograms) posterior PDFs for the galaxy model parameters using constraints 
from the kinematics data set only. The three contours indicate the regions containing, respectively, 68%, 95% and 99.7% of the probability. 



full description of the model also includes a number of additional 
parameters that do not represent physical characteristics of the sys- 
tem (i.e., the line-of-sight inclination, the lens center, the regular- 
ization level, the weights, widths and flattenings of the individual 
MGE Gaussians): these are treated as nuisance parameters and kept 
fixed or marginalized over. 



this case, the joint likelihood can be written simply as 



Pr(^|^,7^)c^exp --^ 



[//Obs _ ,,^od(^0^^^ Nu, Lobs _ ^mod^^^^l^ ^ 

J ^ L J ^ ^^^^ 



1 y [^ii -^2. 



For the likelihood function we follow the standard ap- 
proach of ass uming Gaussian errors on the data points (see e.g. 
iBrewer & Lewis 2006; Suvu et al 2006 : Marshall et aL.2007i\ In 



where the three terms inside the exponential represent the familiar 
misfit functions for the separate contributions of gravitational 
lensing, stellar kinematics and gas kinematics, respectively. We in- 
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Figure 5. Marginalized two-dimensional (contour plots) and one-dimensional (histograms) posterior PDFs for the galaxy model parameters using the combined 
constraints from both lensing and kinematics. The three contours indicate the regions containing, respectively, 68%, 95% and 99.7% of the probability. 



dicate as if^^ the A^^ data points constituting the lensing data set, 
i.e. the pixel values in the galaxy- subtracted observed lensed image 
(see top-right panel of Figure [TJ, each characterized by an uncer- 
tainty (Tij. We denote as i"^^^ the corresponding pixel values of 
the model-predicted image, which are determined by the specific 
choice of model parameters 6 (for example, the bottom left panel 
of Figure [T] shows ^^"^°^(^map), i-e. the model-predicted image in 
the case of the MAP model). An analogous notation (i.e., observed 
values, model-predicted values, uncertainties on the data points) 
applies for the velocity moment 112 in the case of the stellar kine- 
matics and for the circular velocity in the case of gas kinematics. 



The computationally expensive task of exploring and sam- 
pling the joint posterior distribution is accomplished by mak- 
ing us e of the very efficient and robust Mu ltiNest algo- 
rithm (iFeroz & HobsonI l2008l : iFeroz et all l2009h . which imple- 
ments the nested samp ling Monte Carlo technique (ISkillindl200i 
ISivia&Skilli^l2006h . and can provide reliable posterior infer- 
ences even in presence of multi-modal and degenerate multivariate 
distributions. For the analysis of SDSSJ2141, we have launched 
MultiNest with 2000 live points (the live or active points are 
the initial samples, drawn from the prior distribution, from which 
the posterior exploration is started). The large number of live 
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Table 1. Summary of the adopted priors and of the posteriors inferred from the com- 
bined analysis for the model parameters. 



parameter description prior posterior 



^vir/kms ^ 


dark halo virial velocity 


N(255, 45) 


949+44 


r 


dark halo inner logarithmic slope 


U(0, 2) 


0.82:«:g 


C-2 


dark halo concentration 


U(0, 50) 




qh 


dark halo 3D axial ratio 


LAr(l,0.3) 






stellar mass 


U(0, 5) 


1 ^9+0.16 
^•■^^-0.25 


b 


orbital anisotropy parameter: cr^/crj 


U(0, 5) 


1 77+0.30 
' -0.29 



Note. In the prior column, U (a, b) denotes a uniform distribution over the open interval 
{a, b); N(a, b) denotes a normal distribution, with a being the central value and b being 
the standard deviation; LN(a,b) denotes a lognormal distribution, with a being the 
central value for the variable, and b being the standard deviation for the log of the 
variable. In the posterior column, we list, for each parameter, the median value of the 
corresponding marginalized posterior PDF and the uncertainty quantified by taking the 
68% credible interval (i.e., the 16th and 84th percentiles). 



p oints adopted for thi s study (cf., e.g., the MultiNest analysis 
in Barnabe et al ]|20T ll) has allowed us not only to produce very de- 
tailed posterior distributions, but also to gauge the minimum num- 
ber of live points (which is found to be ~200) needed to obtain 
reliable posterior PDFs, which will be very useful in reducing the 
computational load in future analyses of further SWELLS systems. 



7 RESULTS AND DISCUSSION 

In this Section we present and discuss the results of our analysis of 
the disk galaxy SDSS J2141, combining the constraints from both 
the gravitational lensing and the kinematic data sets as described in 
the previous Sections. 

7.1 Inferences on the galaxy model parameters 

As discussed in Section (6] the inferences on the model parameters 
obtained from our analysis are expressed in the form of a multi- 
variate posterior PDF. We consider six free parameters: the virial 
velocity Uyir, inner logarithmic slope y, concentration c_2 and three- 
dimensional axial ratio of the gNFW dark matter halo, the total 
stellar mass M^, and the meridional plane orbital anisotropy ratio b. 
Since visualizing the full six-dimensional surface is challenging, 
we present the inferences in the familiar form of "cornerplots" that 
show all possible one-dimensional and two-dimensional marginal- 
ized posterior PDFs for the six model parameters. The inferences 
obtained when using just one single data set are presented in Fig- 
ures [3] (gravitational lensing only) and |4] (kinematics only), while 
Figure |5] shows the results for the combined lensing and kinematics 
data sets. In each plot, the three contours indicate the regions con- 
taining, respectively, 68%, 95% and 99.7% of the probability, i.e. 
they represent the analogue of the 1,2 and 3cr contours of a Gaus- 
sian distribution. The median value and the corresponding uncer- 
tainties (expressed as 16th and 84th percentiles) for each individual 
parameter are listed in Table [T] 

The constraints provided by kinematics alone are in general 
slightly better than the constraints obtained with a pure gravita- 
tional lensing analysis, in particular for the concentration and the 
stellar mass; obviously, the anisotropy parameter b is completely 
unconstrained in the lensing analysis, and thus, in this case, the pos- 
terior is nothing else than the input uniform prior. The inferences 



on the remaining parameters have uncertainties of similar magni- 
tude in the two cases, but the marginalized posterior PDFs have 
different shapes (note, in particular, the profile for the marginal- 
ized PDF of the halo axial ratio in the two cases), which makes it 
possible to tighten the inferences when lensing and kinematics data 
are considered simultaneously. The effectiveness of the combined 
analysis can be seen in Figure \5\ in particular, we can place tight 
constraints on by clipping both the low-mass and the high-mass 
tails, and we also improve significantly our inferences on the dark 
halo parameters and y, for which we weed out the higher val- 
ues. The meaning and implications of the constraints on the model 
parameters are discussed below in Sections 1712117. 5 1 

The high-probability models drawn from the posterior PDF 
of the combined analysis, and in particular the MAP model, re- 
produce both the lensing and the kinematic observables very accu- 
rately (Figures [T] and 121 respectively). Similarly to what was found 
in Paper II, the most probable lensing models predict a faint coun- 
terimage whose presence is consistent with the noise level. For the 
kinematics, we note that the predicted gas rotation curve matches 
quite well also the data points within the inner 2 arcsec, which were 
conservatively excluded from the fit. 

7.2 Mass budget: baryons and dark matter 

A very intuitive way to visualize the galaxy mass budget as a func- 
tion of radius that is inferred from the combined analysis is pro- 
vided by Figure [6] where we show the circular velocity profile ob- 
tained from the posterior PDF, decomposed into the baryonic and 
dark matter components. The solid lines indicate the median values 
from the posterior PDF, while the shaded regions represent the 68 
per cent confidence intervals. The constraints on the total circular 
velocity v^ot are extremely tight, whereas there are larger uncertain- 
ties on the contributions given by the separate components. Despite 
this, it is clear that the baryonic matter is dominant all over the en- 
tire region for which we have data, with the dark matter component 
becoming progressively more important as we move outwards in 
radius. 

Traditionally, in studies of disk galaxies, the characteristic 
radius at which one measures the dark matter fraction /dm = 
Mdm/Moi is 2.2 times the disk scale length R^, which corresponds 
to the radius at which the circular speed peaks for a razor-thin expo- 
nential disk (see, e.g., .Bershadv et al.,201Q. and references therein). 
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Figure 6. Circular velocity profile inferred from the combined lensing and dynamics analysis. The total circular velocity is shown in black, while the baryonic 
and dark matter components are plotted in blue and in red, respectively. For each component, the solid line represents the median while the shaded region 
encloses 68% of the posterior PDF. The two vertical dashed lines mark, for reference, the location of the Einstein radius (left) and of I.IR^ (right), where 
= 3.58 kpc is the disk scale length of the galaxy. 



In the case of SDSS J2141, we determine /DM(2.2/?d) = 0-28!^ }^, 
by integrating the mass within a spherical radius r = I.IR^. We note 
that this dark matter fraction was inferred from gravitational data 
alone, and is independent of the stellar populations in the galaxy. 
The marginalized posterior PDF for this quantity is shown in Fig- 
ure [T] it is clear that the distribution peaks at around /dm - 0.3; 
dark matter dominated models (i.e. /dm > 0.5), however, are still 
possible, albeit with a low probability of about 9 percent. Models 
without dark matter, on the other hand, are ruled out at more than 
the 3-sigma level, i.e. the probability for /dm < 0.05 is less than 
0.3 per cent. The Paper II analysis of this same galaxy — carried 
out using a less flexible mass model and without including the 
stellar kinematic constraints — found (at lower precision) a higher 
contribution of dark matter at 2.2Rd, i.e. /dm = 0.55^Jj^5, which 
is however still consistent within 1-sig ma with the result that we 
determine here. iTrott et all feOlOh and ISuvu et all (l201lh . by ap- 
plying a combined lensing and dynamics analysis on two differ- 
ent disk galaxies, obtain a fractional amount of dark matter close 
to 45 + 10 percent, which is slightly higher than (but not inconsistent 
with) the /dm that we derive fr om the present analysis. On the other 
hand, van de Yen et al.l feOlOl) . by conducting a combined lensing 
and dynamics stu dy of the same early-type disk galaxy studied by 
iTrott et all tOld) , and adopting a lKroupal (1200 ih IMF, find that the 
upper limit for /dm is only ^ 0.20. Interestingly, the value of /dm 
that we obtain in this analysis is similar to the typical average dark 
matter fraction of about 30 per cent determined for massive early- 
type galaxies within on e effective radius based on lensing and dy- 
namical analysis (e.g., Treu & Koopmans '2 0041 : iTreu et al]|201Ql : 
lAuger et al...2010h : .Spiniello et al.>.2011.) , or by assuming maximal 



stellar component (p-g-, Gerhard et al.ll200lL ICappellari et al.ll200d 



iBarnabe et al ]|20ld[201 1|). However, one should keep in mind that 
/dm has been observed to vary quite significantly between individ- 
ual systems. 

This analysi s also enables us to test whether the "maximum 
disk" hypothesis (Ivan Albada & Sancisi 1986), frequent ly adopted 
in studies of late-type galaxies fe.g. lBell & de Jond200lh , holds for 
the object examine d here. We foll ow the definition of maximum 
disk introduced by ISacketd (Il997h . i.e. Udisk(2.2/?d)/i^tot(2.2/?d) = 
0.85 + 0.10, substituting the circular velocity of the disk Udisk with 
the more relevant circular velocity of the entire baryonic compo- 
nent, Ubar- We find that vi,^(2.2Rd)/vtoi(2.2Rd) = 0.87+[J;g^, which 
corresponds to a maximal disk. From the posterior PDF for this ra- 
tio, the probability that the SDSS J2141 disk is submaximal is about 
10 per cent. 

This result would make SDSSJ2141 something of an out- 
lier when compared with a sample, recen tly studied using dy- 
namical methods, of 30 local disk galaxies jBershadv et al.ll201ll : 
Martinsson 2011). These authors find that, although the ratio 
Udisk(2.2/?d)/i^tot(2.2/?d) increases with the maximum rotation speed 
of the galaxy, even the most massive systems with f;disk(2.2/?d) ^ 
250 km s~^ are submaximal on average. We note, however, that the 
existence of individual massive galaxies consistent with "maxi mal- 
ity" is not ruled out in their study (see, in Figure 2 of Bershady et al.l 
I2OIIL the outlier and the error bars for some of the highest rotation 
velocity systems). In ad dition, it is impor t ant to keep mind that both 
our method and that of iBershadv et al.l (l201lh inevitably rely on 
different assumptions: in our case, for example, a common mass- 
to-light ratio for bulge and disk, in their case assumptions necessary 
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to compare edge-on and face-on galaxies. In addition, the methods 
obtain their information from different parts of the mass distribu- 
tion, with our method being more sensitive to the inner regions, 
owing to the lensing and stellar velocity dispersion constraints. We 
plan to perform a more detailed comparison of the two results once 
data and models for the full SWELLS sample will be available. 

7.3 Constraints on the dark matter halo: shape and profile 



Pure dark matter N-body simulations find that dark halos generally 
have triaxial shapes, with a preference for prolateness, particularly 
in the inner reg ions (e.^ | jing & Sj^ I2OO2I : IXllgood et all l200g : 



lBettetal.112007: 
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Maccio et al.ll2008h~Rece nt numerical work (see, 
]201Q| : IXissera et al.ll201Ql) has shown that includ- 



e.g., 

ing the contribution of the baryons has the effect of modifying the 
overall profile of the dark halo, which flattens to a more axisymmet- 
ric and oblate shape, with an average axial ratio of order 0.85-0.95, 
largely constant with radius. In our study of SDSS J2141, we infer 
from the combined analysis that the dark halo is moderately oblate, 
with an axial ratio = 0.75^^^^. Significantly prolate haloes with 

> 1-25 are strongly disfavoured (i.e., with less than 5 per cent 
probability). This is in good agreement with the numerical results 
on bary on- affected halos, although the median value is slightly flat- 
ter than the typical obtained in the simulations. The axial ratio 
obtained for SDSS J2141 in Paper II (using a less flexible NIE dark 
halo model which does not allow for a variable inner slope), i.e. 

= O.Ql^Q-Jg, was closer to spherical but still consistent, within the 
68% uncertainty, with the more accurate analysis conducted here. 
In contrast, in the only other joint lensing and kinemati cs study of 
a disk galaxy that adopts a non- spherical halo model, ISuvu et al] 
Em ll) find a much more flattened dark matter distribution, with 

= 0.33. These authors adopt a simpler luminous mass model 
than the one considered here (i.e., a razor-thin exponential disk and 
a point-mass bulge) and do not have access to stellar kinematics 
data. 

Including the inner slope y as a free parameter in the dark halo 
mass density model (see Section [3]2]), rather than just adopting a 
fixed isother mal or NEW pro file (as done in pre vious studies, cf. 
Paper II and lSuvu et al.ll201ll : but see also lTrott'et al . 2010, where 
a spherical gNEW halo is used), is important since it allows one to 
account for possible baryon-induced effects, such as adiabatic con- 
traction, that can modify the steepness of the density distribution 
in the galaxy central regions. The data- set at hand, unfortunately, 
does not permit us to place a strong constraint on the inner slope: 
we obtain y = 0.82^^^^, approximately equiprobable over the range 
to 1.5, and perfectly consistent with an unmodified NEW profile. 
The probability that the halo has an inner slope 1.5^7<2isl4 
per cent. We are able to conclude, however, that very steep profiles 
are disfavoured: slopes 1.7 ^ y < 2 have only a 3.5 per cent proba- 
bility, whereas from the adopted uniform prior U(0,2) (see Table[T]) 
one would predict 15 per cent over the same interval. 

We infer a halo concentration parameter c_2 = 7.7^2 5' with 
a low-probability tail for high concentrations (the 95th and 98th 
percentiles fall at c_2 - 17 and c_2 - 30, respectively). The in- 
ferred virial velocity is Uvir = 242^3^. Erom these parameters one 
can derive the posterior PDEs for all other useful quantities char- 
acterizing the gNEW halo, such as the generalized scale radius 
r_2 = 41^^^ kpc, the virial radius Tyir = 315^^3 kpc and the virial 
mass logio(M,ir/Mo) = UASt^.j^. 

Eigure[8] shows a comparison between the dark matter concen- 
tration and the virial velocity from our lensing and dynamics anal- 
ysis of SDSS J2141 (contours), with the predictions from N-body 




Figure 7. Dark matter fraction enclosed within the spherical radius r = 2.2 
disk scale lengths, inferred from the combined lensing and dynamics anal- 
ysis. The median and uncertainty (corresponding to 16th and 84th per- 
centiles) is /dm =0.28+[5 }^. 



simulations in a WMAP 5th year cosmology (iMaccio et al.ll2008r) . 
The uncertainty on our inferred dark matter concentration is quite 
broad ^ 0.2 dex, but is nevertheless in very good agreement with 
the simplest predictions from ACDM (i.e., assuming no contraction 
or expansion of the dark matter in response to galaxy formation). 



7.4 Constraints on the stellar IMF 

The total stellar mass inferred from the combined lensing and dy- 
namics analysis is logjQ(M^/Mo) = 11.12+qq^. This value is very 
well constrained and represents a significant improvement over the 
determination of Paper II, by cutting the low mass tail of the 
posterior PDE of about 0.3 dex. 

In order to draw conclusions on the galaxy IME, we need to 
compare the stellar mass derived from the joint analysis with the 
stellar masses that are inferred from SPS models when assum- 



|Sal 


petei (Il955h 


the 
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the multi-band photometric data set of SDSS J2 141, as described 
in Paper II. However, we note that so far we have neglected the 
contribution of the cold gas: if such a component is present, the 
mass derived above from the combined analysis actually rep- 
resents the total baryonic mass. Therefore, in order to obtain a 
posterior PDE for the stellar mass that can be properly compared 
with the predictions from the SPS models, we need to subtract the 
cold gas fraction, which in disk galaxies (with stellar masses of 
^ IO^^Mq ) accounts for about 20 +10 per cent of the baryonic 
mass (see e.g. lDutton & van den Bosch|[200^ . Under the assump- 
tion that the cold gas is distributed approximately like the stars, 
for each sample in the posterior PDE of we draw a random 
gas fraction /gas e [0, 1] from a Gaussian distribution centered on 
0.2 with a standard deviation of 0.1, and we calculate the quantity 
M^(l - /gas). The gas- subtracted stellar mass derived in this way 
islogio(M^/Mo) = 11.01 



^Q^p about 0.1 dex lower than the value 
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Figure 8. Dark matter concentration c_2 = r^ixlr-i vs dark matter virial 
velocity i^vir- The shaded region shows the prediction (with 1- and 2-sigma 
scatter) from N-bod y simulations in a WMAP 5th year ACDM cosmology 
(iMaccio et aklbOOsh . The contours enclose 68% and 95% of the posterior 
probability from the combined lensing and dynamics analysis. The black 
dot shows the median of the posterior distribution. 



obtained above by ignoring the cold gas contribution. This provides 
a robust lower bound to the stellar mass function. In the future, it 
will be useful to refine these mass estimates by including high res- 
olution constraints on the gas fraction from ALMA observations. 

The posterior PDF for the inferred stellar mass (both with and 
without cold gas) is presented in Figure |9l and compared with the 
distributions obtained for a Chabrier and Salpeter IMF. It is clear, 
just by a visual inspection of this Figure, that our results support 
a Chabrier-like IMF over a Salpeter one. This preference can be 
quantified in a rigorous way by calculating the Bayes factor, i.e. the 
evidence ratio between the two models 



Br 



(26) 



where in our case the likelihood J1{M^) is the posterior PDF for the 
inferred stellar mass, while the priors Pr (M^ | TYc) and Pr (M^ | ) 
are given by the posterior PDFs obtained from SPS models in the 
cases of Chabrier and Salpeter IMFs, respectively. The calculated 
Bayes factor is Bcs = 5.7, which corresponds to substantial evi- 
dence in favour of a Chabrier IM F with respect to a Salpeter IMF 
(see e.g. iKass & Raftervl [l995l , and references therein). In other 
words, if these are the only two possible models, this value of Bcs 
means that there is a 85 per cent probability that the Chabrier model 
is the true one. 

This result corroborates the findin gs of Paper II , and i s 
in agreemen t with the works of, e.g . Bell & de Jond (l200lh , 
iKassin et al.1 (l2006h and Ivan de Yen et al. ( 20l3), which dis- 
favour a Salpeter IMF for disk galaxies, preferring instead 
IMFs t h at pre dict lower stellar masses, such as Chabrier or 
iKroupal (l200lh . Low-mass, fast-rotating early-type galaxies are 
also found to be often inconsistent with a Sa lpeter IMF (e.g. 
ICappellari et al.l2006l : lAuger et al.l20l"0al : lBarnabe et al..201Q) . On 



the other hand, there is m ounting evidence th a t massive ellipti- 
cals favour a Salpeter-like ('Auger e t al. | l2010bl : ItYcu et al.ll20lA 
iBarnabe eTal. 2011: Spiniello et al. 201lh or an even steeper IMF 
( van Dokkum & Conroviboioh . These findings, including the re- 
sults of this work, supp ort the idea tha t the traditional picture of a 
universal IMF (see, e.g., lKroupall2002h might need to be revised in 
favour of a more complicated scenario where the IMF depends on 
the galaxy mass and/or Hubble type. 



7.5 Constraints on the stellar anisotropy 

Determining the shape of the velocity dispersion ellipsoid of disk 
galaxies is important not only in order to understand their global 
dynamical properties, which are related to the formation and evo- 
lution mechanisms, but also because the vertical-to-radial velocity 
dispersion ratio o-^Ictr can be used, together with the galaxy scale 
hei ght distribution , to derive the dynamical mass of the d isk (see 
e.g. lBottemal[l997l : iKregel et al.ll2005l : IWestfall et a0201lh . 

From our analysis, we infer a meridional anisotropy parameter 
b - 1.77^Q29' ^^^^ ^ ^^^y symmetric posterior distribution around 
the median value. In order to facilitate the comparison with the disk 
galaxy studies literature, it is convenient to express the anisotropy 
in the notation - \ - cr\lcr\ (see Sect. 15. lb , where Pz - ^ cor- 
responds to isotropy: in this case we have p^ - 0.43^^^^. These 
results show that, for SDSSJ2141, the velocity dispersion in the 
vertical direction is about three quarters of the radial velocity dis- 
persion. The probability that the velocity dispersion ellipsoid is ap- 
proximately isotropic in the meridional plane (i.e.,/?^ = 0.0 + 0.1) is 
only of order 1 per cent. This confirms that two-integral DF models 
(which are semi-isotropic , i.e. have cr\ = cr\ everywhere, see e.g. 
iBinnev & Tremain3 l2008h do not provide an ideal description of 
the dynamical properties of this galaxy, and a more flexible ap- 
proach allowing for anisotropy, such as the one adopted in this 
work, is warranted. Within the hypothesis of axial symmetry, we 
can then conclude t hat the disk galaxy DF re spects a third integral 
of motion (cf., e.g. jNoordermeer et al.ll2008l) . 

These findings are in agreement with numerous dynamical 
studies of disk galaxies, including the Milky Way, which are well 
known to have velo city dispersion ellipsoids flattene d along the ver- 
tical direction (see Ivan der Kruit & Freemanll201lL and refere nces 
therein). For l ocal d isk galaxies, iGerssen et al. 1 11991 l2000h and 
.Shapiro et all (l2003h determine 0.30 < P^ < 0.75. IWilliams et all 
(l2009h . adopting a dynamical model analogous to the one used 
in this work (i.e., based on anisotropic Jeans equations), find 
0.0 < pj < 0.5 for a sample of 14 spiral and SO galaxies. 
iNoordermeer et al.l (l2008h , using two-dimensional kinematic data 
sets to analyze the dynamics of four early-type disk galaxies, find 
Pz - 0.5, perfectly consistent with the result for SDSS J2141. Re- 
cently, one of the galaxies studied in detail in the Dis ldMass Survey 
was determined to have a more flattened p. = 0.77 dWestfall et all 

Imil). 

We remark that the present study represents the first determi- 
nation of the anisotropy parameter for a disk galaxy well beyond 
the local Universe, at a redshift ziens - 0.14 (a previous combined 
lensing and dynamics study of a disk galaxy at a low er redshift, 
^lens - 0.04, was conducted bv Ivan de Ven et al]|2010l who found 
= 0.1 ± 0.1, consistent with the system being semi-isotropic). 
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Figure 9. Comparison of the stellar mass inferred from the combined lensing and dynamics analysis (black) with the stellar mass determined from photometry 
and stellar population synthesis models, assuming a Chabrier IMF (blue) or a Salpeter IMF (red). The grey shaded histogram shows the posterior PDF for the 
inferred stellar mass when a 20% ± 10% contribution in mass from cold gas in included. In the latter case, the Bayes factor in favor of a Chabrier IMF with 
respect to the Salpeter IMF is 5.7, corresponding to Chabrier being preferred at 85%. 



8 CONCLUSIONS 

We have carried out an in-depth, self-consistent analysis of the 
mass and dynamical structure of the lens disk galaxy SDSS J2141 
at redshift 0.14 by combining the constraints from gravitational 
lensing, Ha rotation curve and stellar kinematics. We have adopted 
a flexible axially symmetric mass model consisting of a gNFW 
dark matter halo and a self-gravitating stellar distribution obtained 
from the MGE parametrization of the observed luminous profile. 
We have modelled the kinematics by means of anisotropic Jeans 
equations in order to allow for a velocity dispersion ellipsoid that 
is flattened in the meridional plane, as is typical for disk galaxies. 

This work improves in several ways (namely, the inclusion of 
stellar kinematics constraints and the upgraded mass and dynami- 
cal model) on the study of this same object described in Paper II, 
and represents the most accurate and detailed analysis to date of 
the dark and luminous mass profile of a disk galaxy beyond the lo- 
cal Universe, i.e. at redshift > 0.1. The main conclusions of this 
analysis can be summarized as follows: 



(i) The spherical dark matter mass fraction within I.IR^ is 
/dm = 0-28^q}q, independent of assumptions on the stellar popu- 
lations in the galaxy. The dark matter fraction increases with ra- 
dius, but does not become dominant within the range probed by the 
observations, which extend to approximately R = \4 kpc. Models 
without dark matter (i.e., /dm < 0.05) are ruled out at more than 
the 3-sigma level. 

(ii) We test the maximum disk hypothesis: we find that, at 2.2Rd, 
the fractional contribution of the baryons to the total circular veloc- 
ity is 0.87^^ i|Q. This corres ponds to a maximal disk (following the 
definition of Isacket3ll997h : the probability of having a submaxi- 
mal disk for SDSSJ2141 is 10 per cent. This is in disagreement 
with recent studies of local disk galaxies (e.g., the DiskMass Sur- 



vey, 'Ber shady et al ] |201lLlMartinssonll201lh . which typically find 
submaximal disks. 

(iii) The gNFW dark matter halo is characterized by a virial ve- 
locity Uvir = 242^39 kms~^ and a concentration parameter c_2 = 
7.7^2 5' implying a generalized scale radius r_2 = 41^^^ kpc. This 
is in very good agreement with the predictions from N-body sim- 
ulations in a ACDM cosmology (i.e., assuming no contraction or 
expansion of the halo in response to galaxy formation). 

(iv) The inner slope of the dark matter halo is only weakly con- 
strained, 7 = 0.82^Q 34, and is consistent with an unmodified NFW 
profile (y = 1). We can still conclude, however, that very steep 
inner profiles with y ^ 1.7 are disfavoured. 

(v) The dark matter halo is moderately oblate, with a three- 
dimensional axial ratio = 0-75^^^^, and a very low probabil- 
ity for significantly prolate haloes (i.e., ^ 1.25). Recent high- 
resolution simulations (e.g. lAbadi et aDl201Ql : iTissera et al.ll201Qh 
find that the baryons have the effect of turning the prolate triaxial 
dark matter halos into roughly oblate spheroids, a scenario that is 
consistent with the results of this work. 

(vi) The total baryonic mass is tightly constrained by the com- 
bined lensing and dynamics analysis, and is determined to be 
logio(M^/Mo) = 11.12!E!:JJ^, independent of the IMF. When ac- 
counting for the expected cold gas contribution, we obtain a stel- 
lar mass logio(M^/Mo) = ILOl^fj-^^. This value is in excellent 
agreement with the stellar mass that is predicted when assuming 
a Chabrier IMF. Model comparison shows that there is substantial 
evidence in favour of a Chabrier IMF with respect to a Salpeter IMF 
(the Bayes factor is 5.7, corresponding to a 85 per cent probability). 

(vii) We infer a meridional anisotropy parameter fi^ - 0.43 ^^^^ 
implying that, for SDSSJ2141, the velocity dispersion ellipsoid 
in the meridional plane is flattened along the vertical direction, in 
agreement with most studies of local disk galaxies. Semi-isotropic 
models (i.e., ^ 0) are ruled out at a very high confidence level, 
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corroborating the evidence that the dynamics of disk galaxies is not 
adequately described by two-integral DFs, and a third integral of 
motion is required. 
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